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QJi [ Abstract 

By analytic continuation to real 9 of data obtained from numerical 
"j^ i simulation at imaginary 9 we study the Haldane conjecture and show 

' that the 0(3) non-linear sigma model with a 9 term in 2 dimensions 

becomes massless at 9 = 3.10(5). A modified cluster algorithm has 
been introduced to simulate the model with imaginary 6. Two differ- 
ent definitions of the topological charge on the lattice have been used; 
■ one of them needs renormalization to match the continuum operator. 

^ , Our work also offers a successful test for numerical methods based on 

Q I analytic continuation. 
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1 Introduction 



Several years ago Zamolodchikov and Zamolodchikov introduced an inte- 

grable S'-matrix for massless particles which was associated to the 2-dimen- 
sional 0(3) non-linear sigma model with a topological 9 term for 6 = n [1]. 

Therefore the 2-dimensional 0(3) sigma model with ^ = tt is possibly 
gapless. Actually this conclusion was previously achieved by Haldane and 
Affleck. They worked out a low energy description of the 1-dimensional chain 
of quantum half-integer spin with antiferromagnctic coupling finding that it 
and the above-mentioned 0(3) sigma model share the same long distance 
properties [2, 3, 4]. Moreover it was found that antiferromagnctic quantum 
spin chains are gapless for half-integer spins [2, 5]. 

Later Affleck and Haldane [6] argued that the critical theory for the half- 
integer quantum antiferromagnctic spin chain is the Wess-Zumino-Witten 
model with a topological coupling k = 1. This model is the stable fixed 
point of the 2-dimensional 0(3) sigma model with a vacuum angle 9 = tt. 

In addition, two numerical calculations of the partition function for the 
0(3) model in the presence of a ^ term [7] yield indications that the the- 
ory undergoes a second order phase transition at ^ = tt (although the two 
analyses disagree about the universality class). 

In the present work we introduce a direct numerical method to verify the 
Haldane conjecture for the 2-dimensional 0(3) non-linear sigma model. The 
idea is to perform a Monte Carlo study on the lattice of the mass gap in the 
model as a function of the 9 parameter and to show that it vanishes at a 
precise value of 9, called 6'cnd) which, following Haldane, should be ^^cnd = 
TT. We overcome the sign problem by simulating the theory at imaginary 9 
and analytically continuing the results to the real 9 values. To this end we 
introduce a new cluster algorithm that works for imaginary non-zero theta. 

2 Lattice implementation 

The continuum expression for the action of the model is 




2 



Q{x) = ^e^''eabcr{x)d,<p\x)d,<P%x) , (1) 
where g is the couphng constant and Q{x) is the topological charge density. 

— * 

(j){x) is a 3-component unit vector that represents the dynamical variable 
at the site x. We have regularized this action on a square lattice by the 
expression 

SL^AL-iOLQi, AL = -(3^$(x)-$(x + fl) , (2) 

where Ql = J2xQl{x) is the total lattice topological charge and Ql{x) is 
the lattice topological charge density. f3 is the inverse bare lattice coupling 
constant and 9l is the bare vacuum angle. In general, 9l ^ and the point 
where the lattice regularized model becomes massless will be called ^L,end- 

The action used in (2) is the simplest one on the lattice. More com- 
plicated actions and expansion parameters boast better scaling and asymp- 
totic scaling properties and hence they are more suited for the calculation of 
masses [8, 9, 10]. However our interest lies only on the vanishing of the mass 
gap at a particular value of and such a property is clearly unaffected by 
the slow convergence of the series. 



3 Choice oi Ql 

Let us discuss the regularization of the topological charge density to be used 
in our Monte Carlo simulation. Q counts how many times the configuration 
of spin variables winds around the unit sphere. Hence Q takes on integer 
values. Configurations with (—1) winding number are called instantons 
(anti-instantons) [11]. 

We have made use of two different lattice regularizations for the topolog- 
ical charge density. The first one [12] 

Q^l\x) = ^e>^''ea,cnx){ct>\x + n)-<l>\x-n))- 

l<l>'{x + <t^\x - V)) , (3) 

is a symmetrical discretization of the expression for Q{x) in Eq.(l). 

The second lattice regularization is defined on triangles (not on single 
sites). Every plaquette of a square lattice can be cut through a diagonal 



3 



into two triangles. If we call (61, (62 and 03 the fields at the sites of the 
three vertices (numbered counter-clockwise) of one of these triangles then 
the fraction of spherical angle subtended by these fields is Q^^\A) which 
satisfies [13] 

exp (27riQf (A)) = ^ (l + 0i • 02 + 02 • 03 + 03 • 0i + ¥1 • (02 X 0";) ) , (4) 
where = 2(1 + 0i ■ (1 + 02 ■ 03) (1 + 03 ■ 0i) and ( A) e [-|, +|]. The 

(2) 

sum of (A) over all triangles yields the so-called geometric topological 
charge Qi\ 

In general, a regularization of Q does not lead to integer values on a 
single configuration. To recover integer results for Ql on ensembles of con- 
figurations that belong to the same topological sector, we must renormalize 
this operator. The lattice and the continuum topological charges are related 
by [14] 

Qr' = zS'''Q, (5) 

(1 2) 

Zq being the corresponding renormalization constant which is UV finite 

since the topological charge operator has no anomalous dimensions in the 

model under study. 

Zq can be calculated either in perturbation theory [14] or by a non- 

perturbative numerical method [15]. We have used the latter. In a nutshell 

it works in the following way: a classical instanton (with topological charge 

+ 1) is put by hand on the lattice and then 100 updating steps are applied 

(we used the Heat-Bath algorithm on the conventional 0(3) non-linear a 

model without a 9 term since the renormalization constant to be used in 

Eq.(5) cannot depend on 6). After every Heat-Bath step the value of Qi''^'' 

is measured and, in order to monitor the background charge and check that it 

(1 2) 

is not varied after the updating step, Q]^' is measured again after 6 cooling 
hits. In the calculation of Zq^ this procedure was repeated 4 ■ 10^ times at 
(3 — 1.5 and 1.6 and 10^ times for (3 — 1.7 and 1.75. The average of Q^'^^ on 
configurations within the topological sector +1 yields Zq''^\ 

The above non-perturbative method is summarized by the expression 



^(^,'^) _ /l-instanton^^ Ql' 6Xp (— 



/l-instanton ^0 CXp (-A^) 
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The restricted path integral runs over all configurations (fiuctuations) that 

(2) 

preserve the background of one instanton. Since the geometric charge Q)^ 
is +1 till the background classical configuration is one instanton (whatever 
the fiuctuations are), the expression (6) yields Z^q'^ = 1 for all P [16]. 

The determination of Zq^ is not so trivial and an example of such an 

evaluation is shown in Fig.l. Measures of Q^^^ on configurations that belong 

to the topological sector +1 attain to a plateau (in general after a few Heat- 
Bath steps) and stay on it for the rest of the updating steps. The height of 



this plateau is the value of Zq\ In Table 1 the results for Zq' at the values 
of P used in the present work are given. 



(1) 
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Figure 1: Data for (Ql^) start at +1 at the 0-th Heat-Bath step and then 
they go down until reaching a plateau. The horizontal line and grey band 
are the value and error respectively of Zq\p = 1.5). 



The relevant consequence of the above considerations for our study is that 

the vacuum angle 6 is related to the corresponding bare parameter by the 



expression 6 — OlZq'"^^ (which implies 6 — 6l when Q^^' is used). 



(2) 
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4 Cluster algorithm for imaginary 



Although the use of the topological charge density Q\^' requires the knowl- 
edge of a renormalization constant, it brings about the advantage that the 
action Sl in (2) can be simulated on the lattice by use of a fast cluster al- 
gorithm. Instead when the geometric charge Q)^ was used, the model was 
updated by a (rather slow) Metropolis algorithm. 

Let us briefly describe the main characteristics of the new cluster algo- 
rithm expressly devised for the present work. The first part of an updating 
step with the usual Wolff algorithm [17] for the standard 0(3) sigma model 
without a 9 term consists in choosing a random unit vector r in such a way 
that every dynamical field can be split in a component parallel to r and 
the rest, (j){x) = {(j){x) ■ rj r + (j)±{x), where (f)±{x) denotes the part of (f){x) 



orthogonal to r. Then the signs of \^(f){x) ■ fj for all x are updated a la 

Swendsen-Wang as in the Ising model [18] . 

— * 

By introducing the above separation for 0(x) in the expression (3) we can 
re-write it as 



+ {f{x + 2)- f) {do,i - do,-i) + (0(x - 2 ) • r) (do,-i - do,i) } , (7) 



where x±l means the site at the position one step forward (backward) in the 
direction "1" starting from site x and the notation dij stands for the 3x3 
determinant (the three components for each vector must be unfold along the 
rows) 



In this fashion the theory at each updating step looks like an Ising model 
in the bosom of an external local magnetic field h{x) because the expression 
in Eq.(7) is linear in (^0 • r^. Recall that all Monte Carlo simulations have 
been performed with an imaginary vacuum angle 9l — +i'&L, (i^l G R)- By 
gathering all contributions of the type shown in Eq.(7) that contain (^(x) • rj 





Ql\x) = ■ f) {di,2 + 0?-l,-2 + d2,-l + rf_2,l) 

+ {${X + 1 ) • f) (cio,-2 - 0?0,2) + (0(X - 1 ) • f) (cio,2 - C?0,-2) 
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at site X one can readily derive the effective magnetic field at this site, 



h{x) = - -—\(p{x) ■f\[di^2 + d-i,-2 + d2-i + d-2,i 

+ '^-1,-1-2 + '^-1+2,-1 + C^l,l+2 + C^l-2,1 

+ C?2,2-l + '^2+1,2 + C^-2,-2+1 + <^-2-l,-2 )■ (9) 

di+k,j (and analogous terms in (9)) are the straightforward generalization of 

the above definition (8) when the site is obtained by shifting two steps from 
the original position x, the first in the direction i and the second in the 
direction k . 

Hence the last step in the updating consists in applying to the above 
expressions an algorithm valid for the Ising model in presence of a mag- 
netic field. In the literature there are two such algorithms, the Lauwers- 
Rittenberg [19] and the Wang [20, 21] methods. After testing their perfo- 
mances and comparing the corresponding decorrelation times with the usual 
(multihit) Metropohs, Heat-Bath and overHeat-Bath, we decided for the 
Wang algorithm. It consists in placing the magnetic field on an extra, fic- 
titious site (called ghost site) that couples to every Ising spin through the 
value of h{x). Using this coupling on the same footing as all other terms in 
the action, the Fortuin-Kasteleyn clusters [22] are arranged by the Hoshen- 
Kopelman algorithm [23] and then updated with the usual | probability. 

Following the proof given in [17], it can be seen that our algorithm also 
satisfies the detailed balance property. 

5 Results 

Operators representing physical states can be built out of an arbitrary num- 
ber of fundamental fields since supposedly the model is not parity invariant 
ioT 9l 0. As the energy gap is given by the mass of a triplet state [24] 
we studied the correlation functions of operators having one 0(3) index as 
quantum number, 

'di{x) =(f{x) , 'd2(x) = $(x) x${x + l). (10) 

Then we calculated the related wall operators by averaging over the Xi co- 
ordinate (as usual L is the lattice size), Wi{x2) = j; J2xi ^i(^) for i = 1, 2. 
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To single out the correct parity mixture for the physical particle and to 
clean the signal from possible excited states, we extracted the triplet mass 
m by using the variational method of Ref. [25] where the mass is obtained 
from the exponential decay of the largest eigenvalue of the correlation matrix 



0.20 



0.15, 
am 
0.1 Oi 



0.05 



-100 -50 _2 50 100 150 

Figure 2: Behavior of the mass gap (in units of the lattice spacing a) as a 
function of Of^. Main plot: circles {(3 = 1.5), up triangles {f3 = 1.6), squares 
{(3 = 1.7) and down triangles {j3 = 1.75) are the data from the simulation at 
imaginary 6l {9\ < 0) by using the Q^l^ lattice regularized topological charge. 
Each continuous line is the result of the extrapolation described in the text 
and the dashed lines enclose the boundary of its error. Inset: the same for 
the Q^'* regularization: squares {j3 = 1.5) and up triangles {j3 — 1.55). In 
this case 9 = 9l and the position of 6* = tt is indicated. 

5.1 Results for Q^}} 

2 • 10^ decorrelated propagators were measured for all values of j3 and 9l. 
In the main plot of Fig. 2 the results for the triplet mass are shown for four 
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values of (3. The extrapolations in this figure were done by using the func- 
tional form (ci + C2 0\)/{1 + Cz0\). We avoided using a functional form 
dictated by some theoretical argument, like one that for real or imaginary 6^ 
goes like m{9L) = Ci(c| — ^|)^/^ which is, up to logarithmic corrections, the 
Renormalization Group prediction, because such an analytic form implicitly 
assumes the vanishing of the mass at a precise value of 6. Instead wc made 
the extrapolations with ratios of polynomials (which are contemporaneously 
both simple and very general functional forms) in order to leave room for any 
behaviour in the vacuum angle. The results of the analytic continuations are 
given in Table 1. The physical value of 6 where the theory becomes gapless is 
given by ^end = ^L,end-^Q • The numbers in the last column are in fair agree- 
ment with the prediction that the model becomes massless when 9 equals n. 
Similar results (and x^) were obtained from degree 2 or 3 polynomials in 9]^, 
while ratios of higher order polynomials proved to be statistically unlikely 
(their was too large). 



Table 1. Values of Zx' and ^end- 



f3 


L 


(^L,cnd)^ 


4^^ 


xVd-o.f. 


n 

"end 


1.5 


120 


111(5) 


0.285(9) 


0.90 


3.00(12) 


1.6 


180 


94(5) 


0.325(6) 


0.45 


3.15(10) 


1.7 


340 


67(3) 


0.380(6) 


1.04 


3.11(9) 


1.75 


470 


56(3) 


0.412(5) 


0.68 


3.08(9) 



The lattice sizes in Table 1 were chosen large enough to meet at = 
the condition L/^ = L ■ am > 10. Once this inequality holds at = 0, it is 
amply realized at the values of 9l where the simulations were performed as 
inferred from Fig. 2. This fact warrants the absence of significant finite size 
effects. 

5.2 Results for Q^l^ 

In this case a Metropolis algorithm was used for updating and 10"'"' indepen- 
dent propagators were measured for each value of 9 (recall that in the present 
case 9l — 9). We report data for two values of (3. They are displayed in the 



9 



inset of Fig.2. The value of 6 where Haldane predicted the closing of the 
mass gap is indicated with an arrow, 6^ = tt^. The numerical results are 
given in Table 2. Comments similar to the case apply to the extrapola- 
tions shown in the figure. Again the results are in fair agreement with the 
conjecture. 

Table 2. ^end for the operator Q^^ 





L 


(^end)^ 


xVd-o.f. 


6'end 


1.5 


110 


10.4(1.0) 


1.72 


3.22(16) 


1.55 


150 


9.7(1.0) 


0.73 


3.11(16) 



By averaging all results for both topological charge operators and assum- 
ing gaussian errors we obtain that the mass gap vanishes at ^end = 3.10(5). 

6 Conclusions 

We have simulated the 0(3) non-linear sigma model in 2 dimensions with an 
imaginary 9 term at several values of the lattice coupling f3. The mass gap 
was measured and extrapolated towards real 9. In all cases the extrapolation 
vanished at a value of 9 compatible with the Haldane conjecture 9 — tt. Our 
result is ^ = 3.10(5) which agrees within errors with the conjecture. This 
value seems very robust as it is independent of the topological charge density 
operator chosen for the simulation. In particular, an operator Q^^ that 
requires a non-trivial renormalization constant leads to the same conclusion 
than another operator (the geometric charge Ql ) ^^^^ does not renormalize. 

A new fast cluster algorithm was purposely introduced to simulate the 
theory with an imaginary 9 term. It works for the operator Ql"*- Instead, 

(2) 

when the geometric topological charge was used, the theory was simu- 
lated by a (rather slow) Metropolis algorithm. 

A salient outcome of our work is the good performance of the analytic con- 
tinuation from imaginary to real 9. No theoretical prejudices were assumed 
in the functional form used in the continuation, apart from the obvious re- 
quirement that it is analytic. This can be justified by a comparison with the 



10 



phase diagram shown in Ref. [26]. In the case of the geometrical charge our 
largest (3 and 9 = 6l (1.55 and 3.5 respectively) are very far from the line 
of phase transitions; as for the Q^^^ case, our largest (3 and 6 — OlZq^ were 
1.75 and 4.1, which again lie very far from any line of singular points. 

The need to perform Monte Carlo simulations at imaginary values of 9 
is actually a blessing in disguise since it forced us to work at very small 
correlation lengths, as can be clearly seen in Fig.2. Had we studied the 
theory directly at real 9 values, we would have met with severe finite size 
effects. 

A key ingredient for the successful extrapolation was to have got data 
from simulations within a wide range of (imaginary) values oi 9l for all j3, 
{^L = -i9L e [0, 10] when was used and = ^ = -i9 e [0, 3.5] 
for Qi , the difference of intervals being due to the effect of the non-trivial 
renormalization that must be applied to the former). All that looks en- 
couraging for the numerical studies based on the analytic continuation with 
respect to a parameter in the theory, such as in QCD with non-zero chemical 
potential [27]. 
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